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Abstract 

We investigate the evolution of the fractions of late type cluster galaxies as a 
function of redshift, using one of the largest, most uniform cluster samples available. 
The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan 
Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created 
using a single automated cluster finding algorithm applied to uniform data from a 
single telescope, with accurate CCD photometry, thus, minimizing selection biases. 
We use four independent methods to analyze the evolution of the late type galaxy 
fraction. Specifically, we select late type galaxies based on: restframe g — r color, u — r 
color, galaxy profile fitting and concentration index. The first criterion corresponds to 



the one used in the classical Butcher-Oemler analyses. The last two criteria are more 
sensitive to the morphological type of the galaxies. In all the four cases, we find an 
increase in the fraction of late type galaxies with increasing redshift, significant at the 
99.9% level. The results confirm that cluster galaxies do change colors with redshift 
(the Butcher-Oemler effect) and, in addition, they change their morphology to later- 
type toward higher redshift — indicating a morphological equivalent of the Butcher- 
Oemler effect. We also find a tendency of richer clusters to have lower fractions 
of late type galaxies. The trend is consistent with a ram pressure stripping model, 
where galaxies in richer clusters are affected by stronger ram pressure due to higher 
temperature of clusters. 
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1. Introduction 

The Butcher-Oemler effect was first reported by Butcher & Oemler (1978, 1984) as an 
increase in the fraction of blue galaxies (/&) toward higher redshift in 33 galaxy clusters over 
the redshift range 0<z<0.54. Butcher and Oemler's work made a strong impact since it showed 
direct evidence for the evolution of cluster galaxies. Much work regarding the nature of these 
blue galaxies followed. Rakos & Schombert (1995) found that the fraction of blue galaxies 
increases from 20% at z=0.4 to 80% at z=0.9, suggesting that the evolution in clusters is even 
stronger than previously thought. Margoniner & De Carvalho (2000) studied 48 clusters in 
the redshift range of 0.03<z<0.38 and detected a strong Butcher-Oemler effect consistent with 
that of Rakos & Schombert (1995). Despite the trend with redshift, almost all previous work 
has reported a wide range of blue fraction values. In particular, in a large sample of 295 Abell 
clusters, Margoniner et al. (2001) not only confirmed the existence of the Butcher-Oemler 
effect, but also found the blue fraction depends on cluster richness. 

Although the detection of the Butcher-Oemler effect has been claimed in various studies, 
there have been some suggestions of strong selection biases in the cluster samples. Newberry, 
Kirshner & Boroson (1988) measured velocity dispersions and surface densities of galaxies in 
clusters and found a marked difference between local clusters and intermediate redshift clusters. 
More recently, Andreon & Ettori (1999) measured X-ray surface brightness profiles, sizes, and 
luminosities of the Butcher-Oemler sample of clusters and concluded that the sample is not 
uniform. The selection bias, thus, could mimic evolutionary effects. Smail et al. (1998) used 
10 X-ray bright clusters in the redshift range of 0.22<z<0.28 and found that the clusters have 
only a small fraction of blue galaxies. The Butcher-Oemler effect was not observed with their 
sample. Similarly, galaxies in radio selected groups are not significantly bluer at higher redshifts 
(Allington-Smith et al. 1993). Garilli et al. (1996) observed 67 Abell and X-ray selected clusters 
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and found no detectable Butcher-Oemler effect until z=0.2. Fairley et al. (2002) studied eight 
X-ray selected galaxy clusters and found no correlation of blue fraction with redshift. Rakos 
& Schombert (1995)'s sample were selected from the catalog compiled by Gunn, Hoessel & 
Oke (1986) using photographic plates taken only in two color bands. The sample thus have a 
possible bias against red, high redshift clusters. In addition to the possible sample selection 
biases, with the exception of Margoniner et al. (2001), the number of clusters in the previous 
works was small, consisting of a few to dozens of clusters. Therefore the statistical uncertainty 
was large. Many authors also noted that cluster-to-cluster variation of the fraction of blue 
galaxies is considerable. The need for a larger, more uniform sample of clusters has been 
evident. 

There have been various attempts to find another physical mechanism causing the large 
scatter which has been seen in almost all previous work. Wang & Ulmer (1997) claimed the 
existence of a correlation between the blue fraction and the ellipticity of the cluster X-ray 
emissions in their sample of clusters at 0.15<z<0.6. Metevier, Romer & Ulmer (2000) showed 
that two clusters with a bimodal X-ray surface brightness profile have an unusually high blue 
fraction value and thus do not follow the typical Butcher-Oemler relation. They claimed that the 
Butcher-Oemler effect is an environmental phenomenon as well as an evolutionary phenomenon. 
Margoniner et al. (2001) found a richness dependence in the sense that richer clusters have 
smaller blue fractions. They claimed that this richness dependence causes a large scatter in the 
blue fraction-redshift diagram. Therefore, it is of extreme interest to explore an origin of the 
scatter in the blue fraction despite the redshift trend. 

At the same time, various studies using morphological information have reported a simi- 
lar evolution effect in cluster galaxies. Dressier et al. (1997) studied 10 clusters at 0.37<z<0.56 
and found a steep increase in the SO fraction toward lower redshift, compared to nearby clus- 
ters studied earlier (Dressier 1980). Couch et al. (1994,1998) studied three clusters at z=0.31 
and found their SO fraction to be consistent with the trend observed by Dressier et al. (1997). 
Fasano et al. (2000) observed nine clusters at intermediate redshifts (0.09<z<0.26) and also 
found an increase in the SO fraction toward lower redshift. It has been proposed that the increase 
in the SO fraction is caused by the transformation of spiral galaxies into SO galaxies through 
a process yet unknown. These studies, however, need to be pursued further, considering that 
most of the previous work was based on morphological galaxy classification by eye. Although 
it is an excellent tool to classify galaxies, manual classification could potentially have unknown 
biases. (A detailed comparison of human classifiers can be found in Lahav et al. 1995). A 
machine based, automated classification would better control biases and would allow a reliable 
determination of the completeness and false positive rate. A further reason to investigate the 
evolution of cluster galaxies is the sample size. The morphological fraction studies of clusters 
reported so far are based on only dozens of clusters. Furthermore, the clusters themselves have 
intrinsic variety in the fraction of blue/spiral galaxies as reported by various Butcher-Oemler 
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and morphological analyses listed above. Since several authors have suggested that the fraction 
of blue galaxies depends on cluster richness, it is important to use a uniform cluster sample, 
preferably selected by an automated method with a well known selection function. 

Various theoretical models have been proposed to explain the Butcher-Oemler effect and 
the increase of the SO fraction. These models include ram pressure stripping of gas (Spitzer & 
Baade 1951, Gunn & Gott 1972, Farouki & Shapiro 1980; Kent 1981, Fujita 1998; Abadi, Moore 
& Bower 1999, Fujita & Nagashima 1999, Quilis, Moore & Bower 2000), galaxy infall (Bothun 
& Dressier 1986, Abraham et al. 1996a, Ellingson et al. 2001), galaxy harassment (Moore 
et al. 1996, 1999), cluster tidal forces (Byrd & Valtonen 1990, Valluri 1993), enhanced star 
formation (Dressier & Gunn 1992), and removal and consumption of the gas (Larson, Tinsley 
& Caldwell 1980, Balogh et al. 2001, Bekki et al. 2002). It is, however, yet unknown exactly 
what processes play major roles in changing the color and morphology of cluster galaxies. To 
derive a clear model explaining the evolution of cluster galaxies, it is important to clarify both 
the Butcher-Oemler effect and the SO increase, using a large and uniform cluster sample in 
conjunction with a machine based morphological classification. 

With the advent of the Sloan Digital Sky Survey (SDSS; York et al. 2000), which is 
an imaging and spectroscopic survey of 10,000 deg 2 of the sky, we now have the opportunity 
to overcome these limitations. The SDSS Cut & Enhance galaxy cluster catalog (Goto et al. 
2002a) provides a large uniform cluster catalog with a well defined selection function. The CCD- 
based, accurate photometry of the SDSS (Fukugita et al. 1996, Hogg et al. 2002, Smith et al. 
2002) and the wide coverage of the SDSS on the sky allow accurate estimation of blue fraction 
with robust local background subtraction. Although the SDSS is a ground based observation, 
the state-of-the-art reduction software and the accuracy of CCD data make it possible to derive 
morphological classification in an automated way (Lupton et al. 2001, 2002). By using the 
SDSS data set, we are able to study one of the largest samples to date — 514 clusters — to 
the depth of M r =-19.44 (h=0.75). 

The purpose of this paper is as follows. We aim to confirm or disprove the existence of 
the Butcher-Oemler effect using one of the largest, most uniform cluster samples. At the same 
time, we hope to shed light on the morphological properties of the Butcher-Oemler galaxies 
using morphological parameters derived from the SDSS data. Finally we investigate the origin 
of the scatter in the galaxy type fraction versus redshift relation, in hope of gaining some 
understanding about the physical processes responsible for the scatter. 

Since several people reported that field galaxies also evolve both morphologically (Schade 
et al. 1996; Brinchmann et al. 1998;Lilly et al. 1998; Kajisawa et al. 2001;Abraham et al. 
2001;) and spectroscopically (Madau et al. 1996; Lilly et al.1996; Hammer et al. 1997; Treyer et 
al. 1998; Cowie et al. 1999; Sullivan et al. 2000; Wilson et al. 2002), it is of extreme importance 
to compare the evolution of cluster galaxies with that of field galaxies to specify a responsible 
physical mechanism. It is possible that the Butcher-Oemler effect and morphological transition 



4 



of cluster galaxies are more commonly happening including the field region of the universe, thus 
a cluster specific mechanism is not responsible for the evolution of galaxies. However, since the 
SDSS spectroscopic data are not deep enough to probe cosmologically interesting time scale, 
we leave it to future work. 

The paper is organized as follows: In Section 2, we describe the SDSS data and the Cut 
& Enhance cluster catalog. In Section 3, we analyze the late type fraction, both spectrally 
and morphologically. In Section 4, we discuss the possible caveats and underlying physical 
processes in the evolution of galaxies. In Section 5, we summarize our work and findings. 
The cosmological parameters adopted throughout this paper are H =75 km s _1 Mpc -1 , and 
(fi m ,fi A ,fi fc )=(0.3,0.7,0.0). 

2. Data 

The galaxy catalog used here is taken from the Sloan Digital Sky Survey (SDSS) Early 
Data Release (see Fukugita et al. 1996, Gunn et al. 1998, Lupton et al. 1999, 2001, 2002, York 
et al. 2000, Hogg et al. 2001, Pier et al. 2002, Stoughton et al. 2002 and Smith et al. 2002 
for more detail of the SDSS data). We use equatorial scan data, a contiguous area of 250 deg 2 
(145.1<RA<236.0, -1.25<DEC<+1.25) and 150 deg 2 (350.5<RA<56.51, -1.25<DEC<+1.25). 
The SDSS imaging survey observed the region to depths of 22.3, 23.3, 23.1, 22.3 and 20.8 in the 
u,g,r,i and z filters, respectively. (See Fukugita et al. 1996 for the SDSS filter system, Hogg et 
al. 2002 and Smith et al. 2002 for its calibration). Since the SDSS photometric system is not 
yet finalized, we refer to the SDSS photometry presented here as u*,g*,r*,i* and z*. We correct 
the data for galactic extinction determined from the maps given by Schlegel, Finkbeiner & 
Davis (1998). We include galaxies to r*=21.5 (petrosian magnitude), which is the star/galaxy 
separation limit (studied in detail by Scranton et al. 2002) in the SDSS data. 

The galaxy cluster catalog used here is a subset of the SDSS Cut & Enhance galaxy 
cluster catalog (Goto et al. 2002a). There are 4638 clusters in the equatorial region (See Kim 
et al. 2002, Annis et al. 2002, Bahcall et al. 2002 and Miller et al. 2002 for other works on 
the SDSS galaxy clusters). Besides the uniformity of the catalog with its well defined selection 
function, the catalog has very good photometric redshifts, 52=0.015 at z <0.3, which enables 
us to use a large sample of clusters (Goto et al. 2002a; see also Gal et al. 2000 and Annis et 
al. 2002 for photometric redshift methods for clusters). We use clusters in the redshift range of 
0.02<z<0.3 and galaxies brighter than M r = —19.44, which corresponds to r*=21.5 at z=0.3. 
Since several authors in previous work claimed that biases in sample cluster selection can mimic 
the evolutionary effect, it is important to control cluster richness well. We use clusters with 
more than 25 member galaxies between M r * = —24.0 and M r = —19.44 within 0.7 Mpc after 
fore/background subtraction, as explained in the next Section. The large areal coverage of the 
SDSS data enables us to subtract the fore/background counts reliably. We can thus control the 
richness of the sample clusters well. The criteria leave 514 clusters in the region. 
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We stress the importance of the uniformity of the cluster catalog. Although the Abell 
cluster catalog (Abell 1958, Abell, Corwin and Olowin 1989) has been used in many studies, 
it was constructed by eye, and is sensitive to projection effects. When comparing clusters 
at different redshifts, it is particularly important to ensure that the data quality and cluster 
selection techniques are comparable, to avoid the introduction of potential selection biases. The 
SDSS Cut & Enhance cluster catalog used here is constructed using only the data taken with 
the SDSS telescope. Also, clusters are detected using a single algorithm (CE) throughout the 
entire redshift range (0.02<z<0.3). Combined with the well controlled richness criteria, our 
cluster sample is not only one of the biggest but also one of the most statistically uniform cluster 
catalogs. To study colors of galaxies in clusters, it is also important to use a cluster catalog 
created without targeting the red sequence of color magnitude relation of cluster galaxies. 
For example, Gladders & Yee(2000) and Annis et al. (2002) use a color filter targeting the red 
sequence of clusters and find galaxy clusters efficiently without suffering from projection effects. 
These techniques, however, can potentially have biases with regards to the colors of detected 
galaxy clusters, since clusters with a strong red sequence is more easily detected. They may 
not, therefore, be ideally suited to a Butcher-Oemler type of analysis. In contrast, the SDSS 
Cut & Enhance cluster catalog does not pick red galaxies selectively, and is therefore more 
suitable for this study. (Note that the the Cut & Enhance method does use generous color 
cuts. Although the color cut is designed to be wide enough to include blue galaxies in clusters, 
therefore it is not completely free from color originated biases.) In previous work, clusters have 
often been detected using data from only one or two color bands. This can introduce a bias, 
since higher redshift clusters are redder and fainter than lower redshift clusters. The SDSS Cut 
& Enhance cluster catalog detects clusters using four bands of the SDSS data (g,r, i and z), 
which minimizes the bias against redshift. 

3. Analysis and Results 

3.1. Fore /Background subtraction 

We compute fractions of late-type galaxies in four different ways. First, we describe the 
statistical fore/background subtraction method we use. In counting galaxies, all galaxies are 
assumed to be at the cluster redshift to calculate absolute magnitudes. Then galaxies only 
between M r * = —24.0 and M r * = —19.44 are used in the analysis. We count the number of 
late- type/total galaxies within 0.7 Mpc from each cluster. Valotto et al. (2001) claimed that 
the global background correction can not correct background contamination appropriately. 
Following the claim, we use a local background correction. The number of late- type/total 
galaxies in fore/background is estimated using an annular area around each cluster with an 
inner radius of 2.1 Mpc and an outer radius of 2.21 Mpc in the same absolute magnitude range. 
The annulus fore/background subtraction enables us to estimate the fore/background locally, 
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minimizing variations in galaxy number counts due to the large scale structure. When an outer 
annulus touches the boundary of the region, a fore/background count is globally subtracted 
using galaxy number counts in the entire 400 deg 2 region by adjusting it to the angular area 
each cluster subtends. This fore/background subtraction is used in the analyses described in 
sections 3.2-3.5. The fraction of blue/late- type galaxies and its error are computed according 
to the following equations. 

ATlate ATlate 

f _ ly c+f iV / m 
J late njall ATall ' V 1 / 



N late + N late N f f + N f 2(y/ N l f e X Nf + y£ ' Njf } X N% f ) 

1 ' (Nffi - Nf te ) 2 + (Ng f - Nf ) 2 (Nffi - Nf te ) (Ng f - Nf ) { } 

where Nf^f and Nf e represent numbers of blue/late- type galaxies in a cluster region and a 
field region, respectively. N%£f and Nf represent numbers of all galaxies in a cluster region and 
a field region, respectively. The equation (2) assumes that N late and N aU are not independent. 
We explain the derivation of equation (2) in Appendix 2. 

3.2. Blue Fraction 

The blue fraction of galaxy clusters (fb) is measured as the fraction of galaxies bluer in 
g — r rest frame color than the color of the ridge line of the cluster by 0.2 mag. This color 
criterion is equivalent to Butcher & Oemler's (1984) 0.2 mag in B — V and Margoniner et al.'s 
(2000, 2001) 0.2 mag in g — r. The color of the ridge line is measured from the color-magnitude 
diagram using the same color-magnitude box used in measuring photometric redshift (Goto et 
al. 2002a). The colors of the ridge lines are confirmed to agree with empirical color of elliptical 
galaxies observed in the SDSS at the same redshift with less than 0.05 difference in g — r color 
(Eisenstein, private communication). We also use Fukugita et al.'s (1995) model of an elliptical 
galaxy and a galaxy bluer than it by 0.2 mag in g — r. By redshifting these two galaxies, we 
measure S(g — r) in the observed frame, which corresponds to the restframe S(g — r)=0.2 . In 
calculating f , we count galaxies within 0.7 Mpc from the center of each cluster, which is the 
same radius as Margoniner et al. (2000, 2001), and corresponds to the average radius of Butcher 
& Oemler (1984). (We explore possible caveats in using fixed radius in Appendix.) Galaxies 
between M r * = —24.0 and M r * = —19.44 are counted, which corresponds to r*=21.5 at 2=0.3 
for an average /c-correction of all types of galaxies (Fukugita et al. 1995). Compared to the 
field luminosity function of Blanton et al. (2001), this includes galaxies as faint as M*+1.36. 
Fore/background galaxies are statistically subtracted in the way described in Section 3.1. The 
lower left panel of Figure 1 shows /& as a function of redshift. The error in /& is estimated using 
equation (2) and the median values of the errors in fb and z are shown in the upper left corner 
of the plot. Dashed line shows the weighted least-squares fit to the data. Solid lines and stars 
show the median values of the data. The scatter is considerable, but both of the lines show a 
clear increase of fb toward higher redshift. The Spearman's correlation coefficient is 0.238 with 
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significance of more than 99.99% as shown in Table 1. The correlation is weak, but of high 
significance. The lower left panel of Figure 2 further clarifies the evolution effect. A dashed 
and a solid line show normalized distributions of fb for clusters with z<0.15 and 0.15<z<0.3, 
respectively. The two distributions are significantly different at the 98% level, as determined 
by a Kolomogorov-Smirnov test. The slope shown with the dashed line rises up to fb ~=0.2 at 
z=0.3 (look back time of ~3.5 Gyr), which is consistent with previous work such as Butcher 
& Oemler (1978,1984), Rakos & Schombert (1995), Margoniner et al. (2000, 2001), within the 
scatter. We conclude that the Butcher-Oemler effect is seen in the SDSS Cut & Enhance galaxy 
cluster catalog. 

We caution readers on the systematic uncertainties in measuring Marzke et 

al.(1994,1997,1998), Lin et al. (1999) and Blanton et al. (2001) showed that luminosity func- 
tions of galaxy clusters depend significantly on galaxy type, in such a way that the bright end 
of the cluster luminosity function is dominated by redder galaxies and the faint end is dom- 
inated by bluer galaxies. Boyce et al.(2001) and Goto et al. (2002b) showed that a similar 
tendency exists for cluster luminosity functions. This difference in luminosity functions leads 
to a different blue fraction depending on the absolute magnitude range used. Furthermore, if 
the radial distributions of blue and red galaxies are different (e.g. Kodama et al. 2001), the 
fb measurement depends heavily on the radius. When comparing with previous work, there- 
fore, it is important to take account of the exact method used to calculate We discuss the 
uncertainty in measuring blue fractions further in Section 4. 

3.3. Late Type Fraction Using u — r <2.2 

Recently Shimasaku et al. (2001) and Strateva et al. (2001) showed that the SDSS u — r 
color correlates well with galaxy morphologies. In this section we use u — r color to separate 
early(n — r >2.2) and late (u — r <2.2) type galaxies as proposed by Strateva et al. (2001). 
Note that although u — r color is claimed to correlate well with galaxy types, it is still a color 
classifier and thus different from the morphological parameters we investigate in the following 
two sections. The methodology used to measure late type fraction is similar to the one we use 
to measure fb- We regard every galaxy with u — r <2.2 as a late type galaxy. We define f u - r 
as the ratio of the number of late type galaxies to the total number of galaxies within 0.7 Mpc. 
Fore/background subtraction is performed in a way described in Section 3.1. The upper left 
panel of figure 1 shows / n _ r as a function of redshift. The error in f u _ r is calculated using 
equation (2) and the median values of the errors in f u _ r and z are shown in the upper left 
corner of the panel. The dashed line shows the least square fit to the data. The solid lines and 
stars show the median values of the data. As in the case of fb, the scatter is considerable, but 
the weak increase of the late type galaxies is seen. The Spearman's correlation coefficient is 
0.234 and is inconsistent with zero at greater than 99.99% confidence level (Table 1). Again, 
weak but significant correlation is found. The upper left panel of Figure 2 shows distributions 
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of f u - r for z<0.15 clusters and 0.15<z<0.3 clusters with a dashed and solid line, respectively. 
A Kolomogorov-Smirnov test shows that the distributions are different with more than a 99% 
significance. In addition to the increase in /& shown in the last Section, the increase in f u _ r 
provides further evidence of color evolution of cluster galaxies. Furthermore, since u — r color 
of galaxies is sensitive to a galaxy's morphology as shown in Figure 6 of Strateva et al (2001), 
it suggests possible evolution of morphological types of galaxies as well. We investigate the 
morphological evolution of galaxies in clusters in the next subsection. 

3.4- Late Type Fraction Using Profile fitting 

One of the purposes of this paper is to determine if there is a morphological change of 
galaxies in clusters as a function of redshift. The SDSS photometric pipeline (PHOTO; Lupton 
et al. 2002) fits a de Vaucouleur profile and an exponential profile to every object detected in 
the SDSS imaging data and returns the likelihood of the fit. By comparing the likelihoods of 
having an exponential profile against that of a de Vaucouleur profile, we can classify galaxies 
into late and early types. In this section, we regard every galaxy that has an exponential 
likelihood higher than a de Vaucouleur likelihood in r band as a late type galaxy. A galaxy 
with higher de Vaucouleur likelihood in r band is regarded as an early type galaxy. We define 
f exp in the same way as in previous subsections, i.e. f exp is the ratio of the number of late 
type galaxies to the total number of galaxies within 0.7 Mpc. Fore/background counts are 
corrected using the method described in Section 3.1. The resulting f exp is plotted in the lower 
right panel of Figure 1. The error in f exp is estimated using equation (2) and the median values 
of the errors in f exp and z are shown in the upper left corner of the plot. The dashed line 
shows the weighted least-squares fit. The solid lines and stars show the median values of the 
data. The scatter is considerable, but we see the increase of f exp toward the higher redshift. 
The Spearman's correlation coefficient is 0.194, which is inconsistent with zero at more than a 
99.99% confidence level (Table 1). The upper right panel of Figure 2 shows the distributions of 
clusters with z<0.15 and with 0.15<z<0.3 with a dashed and solid line, respectively. The two 
distributions show a difference of more than 99% significance in a Kolomogorov-Smirnov test. 
We emphasize that the galaxy classification used here is purely morphological — independent of 
colors of galaxies. The fact that we still see the increase of the late type galaxies toward higher 
redshift suggests that these Butcher-Oemler type blue galaxies also change their morphological 
appearance as well as their colors. We also point out that the slope of the change is similar to 
the panel in the lower left of Figure 1, which is ~30% between z=0.02 and z=0.3. We note that 
there is a potential bias associated with the use of r band profile fitting throughout the redshift 
range, since the r band wavelength range at z=0.3 is almost that of g band at restframe. We 
investigate this effect in Section 4, and conclude that it is small. Like the blue fraction, the 
morphological late-type fraction is also sensitive to the magnitude range considered. Binggeli 
et al. (1988), Loveday et al. (1992), Yagi et al. (2002a,b) and Goto et al. (2002b) reported 
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luminosity functions of elliptical galaxies have brighter characteristic magnitudes and flatter 
faint end tails compared to those of spiral galaxies in both field and cluster regions. Careful 
attention should be paid to the magnitude range used in an analysis when fractions of spiral 
galaxies are compared. We discuss the uncertainty further in Section 4. 

3.5. Late Type Fraction Using Concentration Parameter 

As another morphological galaxy classification method, we use the inverse of the con- 
centration parameter (Cj„) advocated by Shimasaku et al. (2001) and Strateva et al. (2001). 
We define Cj n as the ratio of Petrosian 50% radius to Petrosian 90% radius in r band. They 
are the radii which contain 50% and 90% of Petrosian flux, respectively. (See Stoughton et 
al. 2002 for more details of Petrosian parameters). Since Ci n is the inverse of a conventional 
concentration parameter, spiral galaxies have a higher value of Ci n . Following Strateva et al. 
(2001), we use Cj n =0.4 to divide galaxies into early and late type galaxies. Readers are referred 
to Morgan (1958,1959), Doi, Fukugita & Okamura (1993) and Abraham et al. (1994, 1996) 
for previous usage of concentration of light as a classification parameter. fa n is defined as the 
ratio of the number of galaxies with C, n >0.4 to the total number of galaxies within 0.7 Mpc as 
in the previous subsections. Note that our early type galaxies with C in <0.4 include SO galaxies 
in addition to elliptical galaxies since discerning elliptical and SO galaxies is very difficult with 
the SDSS data, in which the seeing is typically 1.5". (See Shimasaku et al. 2001 and Strateva 
et al. 2001 for the correlation of Cj„ with an eye classified morphology). Fore/background 
number counts are corrected as described in Section 3.1. The absolute magnitude range used 
is —24 < M r < —19.44. The upper right panel of Figure 1 shows fa n as a function of redshift. 
Since the classification using Cj n =0.4 leans toward late type galaxies, the overall fraction is 
higher than the other panels in the figure. The increase of late type fraction, however, is clearly 
seen. The dashed line shows the weighted least-squares fit. The solid lines and stars show the 
median values of the data. The error in fcin is estimated using equation (2) and the median 
values of the errors in fc in and z are shown in the upper left corner of the plot. The Spearman's 
correlation coefficient is 0.223 with significance of more than 99.99% as shown in Table 1. The 
upper right panel of Figure 2 further clarifies the evolution effect. The distribution of z<0.15 
clusters in a dashed line and the distribution of 0.15<z<0.3 clusters in a solid line show a 
difference with more than a 99% significance level. We stress that the galaxy classification 
based on this concentration parameter is purely a morphological one. In this morphological 
classification, we still see the increase of the late type fraction just like the increase of fb — as 
if observing the morphological equivalence of the Butcher-Oemler effect. The increase in fcin 
combined with the increase in f exp provides rather firm evidence of morphological change in 
the Butcher-Oemler type galaxies. Possible caveats in the usage of Cj n and comparisons with 
previous works are discussed in Section 4. 



10 



3. 6. On the Origin of the Scatter 

In the last four sections, we observed the increase of late type fractions toward higher 
redshift in all the four cases. At the same time, we see a significant amount of scatter around 
the late type fraction v.s. redshift relations. Although the errors on these measurements are 
also large, this scatter might suggest that there might be one or more physical properties which 
determine the amount of late-type galaxies in clusters. Table 3 compares median error sizes 
°f hi fu-r, fexp and fcin with scatters around the best-fit lines (dotted lines in Figure 1). In 
fact, in all cases, real scatters are larger than the statistical errors. In the literature, several 
correlations are proposed such as with X-ray shapes of clusters (Wang et al. 2001, Metevier et al. 
2000), and cluster richness (Margoniner et al. 2000). Our cluster richnesses are plotted against 
redshift in Figure 3. Richnesses are measured as numbers of galaxies between M r = —24.0 and 
M r = —19.44 within 0.7 Mpc after fore/background subtraction, as explained in Section 3.1. In 
figure 3, this richness has no apparent bias with redshift. In Figure 4, the difference of the late 
type fraction from the best-fit line is plotted against cluster richness. The circles and solid lines 
show median values. In all the panels, there is a clear tendency of richer clusters having a lower 
fraction of late type galaxies. This tendency is in agreement with Margoniner et al. (2000) 
who found richer clusters had lower blue fractions. We further discuss the richness dependence 
of the late type fraction in Section 4 and Appendix. As an alternative parameter to X-ray 
shape, we plot the difference from the best-fit line against cluster elongation in Figure 5. The 
elongation parameter is taken from Goto et al. (2002a), which measured the ratio of major and 
minor axes in their enhanced map to find clusters. Circles and solid lines show median values. 
No obvious trend is seen here. Our result seems to agree with Smail et al.'s (1997) caution that 
a correlation between ft, and cluster ellipticity found by Wang et al. (1997) could be due to a 
small and diverse sample. However, since distribution of galaxy positions might not represent 
cluster ellipticities well, we do not conclude that there is no dependence on cluster ellipticities. 
The dependence should be pursued further in the future, ideally using X-ray profile shape with 
a large sample of clusters. 

4. Discussion 

4-1. Morphological k- correction 

In the upper right panel of Figure 1, we use Cj„ (inverse of concentration index) in the 
r band to classify galaxies throughout our redshift range (0.02<z<0.3). This could potentially 
cause redshift dependent biases in our calculation of Ci n . Since the universe is expanding, 
by analyzing the observed r band data, we are analyzing bluer restframe wavelengths in the 
higher redshift galaxies. In fact, the r band at z=0.3 is almost g band in the restframe. 
Various authors have pointed out that galaxy morphology significantly changes according to 
the wavelength used (e.g. Abraham et al. 2001). To estimate how large this bias is, we plot the 
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normalized distributions of Cj n in g and r bands in Figure 6 using the galaxies with 0.02<z<0.03 
in the SDSS spectroscopic data (1336 galaxies in total; See Eisenstein et al. 2001, Strauss et 
al. 2002 and Blanton et al. 2002 for the SDSS spectroscopic data). In this redshift range, the 
color shift due to the expansion of the universe is small. We use this to study the dependence 
of the Ci n parameter on the restframe wavelength. At z—0.3, r band corresponds to restframe 
g band. The solid and dashed lines show the distribution for g and r bands, respectively. The 
two distributions are not exactly the same, but the difference between the two distributions is 
small. We summarize the statistics in Table 4. There are 802/1336 galaxies with Ci n >0.4 in 
g band, and 787/1336 galaxies have Cj„ > 0.4 in r band. The difference is 15/1336 galaxies, 
which is 1.1% of the sample. In section 3.5, the change in fa n is ~30%. The effect of the 
morphological k-correction is therefore much smaller. We also point out that this analysis 
assumes the largest difference in redshift (0.02<z<0.3), therefore it gives the upper limit of the 
bias. Since the majority of our clusters are at z ~0.2, the wavelength difference between the 
observed and restframe bands is typically much smaller. We conclude that the effect of the 
morphological k-correction is much smaller than the change in /p m we observed in Section 3.5. 

In Section 3.4, we use the r band fit for all galaxies in our sample. The same redshift effect 
could potentially bring bias to our analysis. In Table 5, we limit our galaxies to 0.02<z<0.03 
and count the fraction of late type galaxies in the g and r bands corresponding to observed 
frame r band at z=0.0 and z=0.3. We list the number of galaxies with exponential likelihood 
higher than de Vaucouleur likelihood in column 1, the total number of galaxies in column 2, 
and the ratio of columns 1 to 2 in column 3. As shown in the 3rd row, the difference in the 
fraction of late type galaxies between g band data and r band data is only 2.5%, which is much 
smaller than the f exp change we see in the upper right panel of Figure 1 (~30%). We conclude 
that the change of f exp and fcin is not caused by the small redshift bias in using r band data 
throughout the redshift range. 

4-2. Seeing Dependence 

Another possible source of bias in measuring f exp and fa n is the dependence on the 
seeing, relative to the size of the galaxies. At higher redshift, the size of a galaxy is smaller and 
a seeing convolution could be more problematic. Especially for the concentration parameter 
(Ci n ), galaxy light becomes less concentrated when the seeing size is comparable to the galaxy 
size, and thus, the effect could cause a bias towards higher C in values. To check this, we plot 
fan against the point-spread function (PSF) size in the r band for two redshift limited samples 
in Figure 7. Open squares and solid lines show the distribution and medians of low z clusters 
(z<0.15). Filled triangles and dashed lines show the distributions and medians of high z clusters 
(0.24<z<0.3). For the median measurements, bins are chosen so that equal numbers of galaxies 
are included in each bin. 1 a errors are shown as vertical bars. As expected, lower redshift 
clusters show almost negligible dependence on seeing size. Higher redshift clusters show about 
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a 5% increase in fa n between the best and worst seeing size. The evolution effect we see in the 
upper right panel of Figure 1 is more than 20%. Furthermore, as the distribution of seeing in 
Figure 8 shows, 87% of our sample galaxies have seeing better than 2.0" . Therefore we conclude 
that varying seeing causes a small bias which is significantly weaker than the evolution we find 
in Section 3. The effect of varying seeing is less significant for the f exp parameter. In Figure 9, 
we plot f exp against seeing size for two redshift samples with redshift ranges and symbols as in 
Figure 7.1a errors, shown as vertical bars, are dominant. There is no significant correlation 
of fexp with seeing size. 

4-3. Radius, Fore/background Subtraction and Cluster Centroids 

Throughout the analyses in Section 3, we use a 2.1-2.21 Mpc annular region for 
fore/background subtraction. In return for taking cosmic variance into account, annular (lo- 
cal) background subtraction has larger errors than global background subtraction due to its 
smaller angular area coverage. However the difference is not so large. In case of blue galaxy 
counts (/&) in the background, the median Poisson (1 a) uncertainty for global background is 
12.2%, whereas 1 a variation of local background is 12.6%. This increases the errors, but only 
by 0.4 points. The actual effect to the late-type fraction is plotted in Figure 10. Solid lines 
show distributions for our default choice of 0.7 Mpc radii and 2.1-2.21 Mpc annuli. Dashed 
lines show distributions for global fore/background subtraction, where fore/background sub- 
traction is performed using global number counts of galaxies for all the clusters in the sample. 
A Kolomogorov-Smirnov test between two samples does not show any significant difference. 

For cluster radius, we use 0.7 Mpc, since we do not have information about the virial 
radii of each system. It is, however, ideal to use virial radii since, for example, in a standard 
cold dark matter cosmology, virial radii at a fixed mass scales as cx (1+z)" 1 . Another pos- 
sible cause of uncertainty is the accuracy in deciding cluster centers. In this work, a center 
position of each cluster is taken from Goto et al. (2002a), and is estimated from the position 
of the peak in their enhanced density map. Although, from Monte-Carlo simulations, cluster 
centroids are expected to be determined with an accuracy better than ~40 arcsec, the offsets 
have a possibility to introduce a bias in our analyses. We test different choices of these param- 
eters in Figure 10. Dotted lines show distributions where radii change as 0.7 x (l+z)^ 1 Mpc 
assuming a standard cold dark matter cosmology. Long dashed lines show distributions when 
the position of brightest cluster galaxy (within 0.7 Mpc and Mr <-24.0) is used as a cluster 
center. Kolomogorov-Smirnov tests show no significant difference in any of the above cases. 
In all cases, the probability that the distributions are different is less than 26%. Our results 
in Section 3 are thus not particularly sensitive to our choice of annuli, radii or cluster centers. 
We further pursue the effect of radius dependence of blue/late type fractions in Appendix, and 
show that it does not change our main results. 
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4-4- Comparison with Late-type Fraction from Spectroscopy 

To further test our late-type fraction measurement, we compare the late- type fraction 
obtained from the SDSS spectroscopic data with that obtained from the SDSS imaging data 
in figure 11. Since the SDSS spectroscopic data are limited to r* < 17.77 (Strauss et al. 2002), 
the comparison can be done only for clusters with z<0.06. In the literature, three clusters are 
found to satisfy these criteria in the region used in this study. The clusters include ABELL 
295, RXC J01 14.9+0024, and ABELL 957. For these clusters, late- type fractions are measured 
in the same way as in section 3. Late-type fractions from spectroscopy are measured using 
galaxies within 0.7 Mpc and <5z=±0.005 from the redshift of each cluster. Note that there 
is no fore/background correction for spectroscopic late- type fraction. In figure 11, all points 
agree with each other within the error. The good agreement suggests that our fore/background 
subtraction technique described in section 3.1 works properly. It would be ideal to perform the 
same test for high redshift clusters as well. However, the SDSS spectroscopic data are not deep 
enough to perform the test for higher redshift clusters. 

4-5. The Butcher- Oemler Effect 

The Butcher-Oemler effect — an increase in the ratio of blue galaxies in clusters as a 
function of redshift — is strong evidence of direct evolution of the stellar populations in galaxies; 
it has been studied by numerous authors in the past. In this section, we compare our results with 
previous work. Since different authors use different cluster samples, color bands, cosmology, 
absolute magnitude ranges and methods of fore/background subtraction, which could affect the 
comparison, we emphasize the differences in analysis by each author. Note that one important 
difference is that some previous work used a sample of quite rich clusters, e.g. clusters with 
more than 100 members in magnitude and radius ranges comparable to those adopted in this 
study. Poorest systems in our sample have only 25 member galaxies after fore/background 
subtraction. Thus, difference in cluster samples could cause a difference in results. 

Butcher & Oemler (1978, 1984) studied 33 clusters between z=0.003 and z=0.54. They 
used galaxies brighter than My=-20 (ft.=0.5 and <7o=0.1) within the circular area containing 
the inner 30% of the total cluster population. They found ft, increases with redshift for z > 
0.1. Their /& at z=0.3 is ~0.15, which is slightly lower than our value. Considering the large 
scatter in both their and our samples, we do not claim that our results are inconsistent with 
their value. Note that Andreon & Ettori (1999) found a trend of increasing X-ray luminosity 
with increasing redshift in the sample clusters of Butcher & Oemler (1984). 

Rakos & Schombert (1995) studied 17 clusters using Stromgren uvby filters. Due to the 
usage of the narrow band filters redshifted to the cluster distance, their study did need to use 
model-dependent k- correct ions. However, their high-redshift cluster sample is drawn from that 
of Gunn, Hoessel & Oke (1986) which is based on Illa-J and IIIa-F photographic plates. At 
z > 0.5, these plates measure the rest-frame ultraviolet to blue region of the spectrum. Thus 
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the cluster catalog will be biased toward clusters rich in blue galaxies at high redshift. Rakos 
&: Schombert found /& ~0.25 at z=0.3, which is slightly higher than the estimation of Butcher 
& Oemler (1984) but in agreement with our results. 

Margoniner et al. (2000) studied 44 Abell clusters between z=0.03 and z=0.38. They 
used galaxies between M r = —21.91 and M r = —17.91 (h = 0.75) within 0.7 Mpc of the clus- 
ter center. The fore/background counts are subtracted using five control fields. Their re- 
sults are more consistent with the steeper relation estimated in 1995 by Rakos and Schombert 
than with the original one by Butcher and Oemler in 1984. The results are also consistent 
with ours. Margoniner et al. (2001) extended their work to 295 Abell clusters and found 
/6=(1.34±0.11) xz-0.03 with a r.m.s. of 0.07, which is in agreement with our fitted function 
shown in Figure 1. 

Ellingson et al. (2001) studied 15 CNOC1 clusters (Yee, Ellingson, & Carlberg 1996) 
between z=0.18 and z=0.55. Since they used spectroscopically observed galaxies, they do not 
suffer from the fore/background correction (but see Diaferio et al. 2001). They used galaxies 
brighter than M r = —19.0 within r2oo from the cluster center (with an average of 1.17/t -1 Mpc). 
Their best fit shows fb ~0.15 at z=0.3. The scatter in their Fig. 1 and our data are both 
substantial. Thus, we can not conclude this is inconsistent with our results. 

All of these authors found considerable scatter in ft v.s. z plot as is seen in our Figure 
1. It is promising that our results are consistent with the previous work within the scatter, 
despite the differences in the radial coverage and magnitude ranges used. 

4-6. Morphological Butcher- Oemler effect 

In sections 3.3, 3.4, and 3.5, we found an increase in the fraction of late type galaxies 
selected by morphological parameters with increasing redshift — as if the Butcher-Oemler effect 
is happening morphologically. Perhaps revealing this morphological Butcher-Oemler effect is 
the most striking result of this study. It suggests that the Butcher-Oemler blue galaxies change 
their morphology from late to early type at the same time that they change their color from 
blue to red. Although accurately quantifying the fraction of galaxies which experience the 
morphological Butcher-Oemler effect is difficult due to the considerable scatter in the data, our 
best-fit lines suggest that ~30% of galaxies in clusters go through this transition between z=0.3 
and z=0.02. 

In previous work, Dressier et al. (1997) found a deficit of SO galaxies in 10 intermediate 
(z~0.5) clusters by classifying galaxy morphology in the HST image by eye. They claimed that 
many SOs needed to be added to reach the fraction of SOs found in present clusters (Dressier 
1980). Couch et al. (1994, 1998) also found an indication of morphological transformation 
in the Butcher-Oemler galaxies by studying three rich clusters at z=0.31. Later, Fasano et 
al. (2000) showed that spiral galaxies are, in fact, turning into SO galaxies by observing nine 
clusters at intermediate redshifts and analyzing them together with higher redshift clusters in 
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the literature. Their galaxy morphology was also based on eye classification. Our SDSS data is 
taken using ground based telescopes with moderate seeing (~1.5" ), and thus does not allow us to 
separate SO galaxies from elliptical galaxies as the HST does. The advantage of our classification 
is its automated nature, which allows accurate reproducibility and quantification of systematic 
biases. In particular, it is easy to compute the completeness and contamination rate for the 
automated classification, based on simulations; for the present sample, the completeness and 
contamination rate of the parameters are given in Shimasaku et al. (2001) and Strateva et 
al. (2001). Furthermore, an automated galaxy classification is easier to reproduce in future 
observational work and in detailed computer simulations. Although we can not distinguish SOs 
from ellipticals, the increase of blue fraction and increase of late type galaxies toward higher 
redshift is qualitatively consistent with the process of SO production over the interval in cosmic 
time suggested by previous investigations. 

Various physical mechanisms could be the cause of the morphological and spectral 
Butcher-Oemler effects. Possible causes include ram pressure stripping of gas (Gunn & Gott 
1972, Farouki & Shapiro 1980; Kent 1981, Abadi, Moore & Bower 1999, Quilis, Moore & Bower 
2000), galaxy infall (Bothun & Dressier 1986, Abraham et al. 1996a, Ellingson et al. 2001), 
galaxy harassment via high speed impulsive encounters (Moore et al. 1996, 1999), cluster tidal 
forces (Byrd & Valtonen 1990, Valluri 1993) which distort galaxies as they come close to the 
center, interaction/merging of galaxies (Icke 1985, Lavery & Henry 1988, Bekki 1998), and 
removal & consumption of the gas due to the cluster environment (Larson, Tinsley & Caldwell 
1980, Balogh et al. 2001, Bekki et al. 2002). Mamon (1992) and Makino & Hut (1997) showed 
that interactions/mergers can occur in a rich cluster environment despite the high relative ve- 
locities. Shioya et al. (2002) showed that the truncation of star formation can explain the 
decrease of SO with increasing redshift. It has been known that preheating of intergalactic 
medium can effect morphologies of galaxies by strangling the gas accretion (Mo & Mao 2002; 
Oh & Benson 2002). In fact, Finoguenov et al. (2003) found the filamentary gas in Coma 
cluster and predicted quiescent star formation in galaxy disks around the filament. Although 
our results provide some important clues, pinpointing what processes are responsible in the 
morphological and spectral Butcher-Oemler effect is a more difficult challenge. 

Our results suggest that the cause will be one that affects both color and morphological 
appearance of galaxies at the same time. Couch et al. (1998), Dressier et al. (1999) and 
Poggianti et al. (1999) found "passive spirals", which are galaxies with spiral morphology but 
without star formation. They are probably belong to the same population as "anemic spirals" 
found by van den Bergh (1976). The mechanism creating "passive spirals" or "anemic spirals", 
however, affects only the color of galaxies and, thus, probably is not the main mechanism that 
accounts for the entire effect. The increase of morphologically late type galaxies toward higher 
redshifts at the same time as the increase of blue galaxies is consistent with mechanisms which 
affect the gas supply (e.g. ram-pressure stripping, galaxy infall). However if the infalling rate 
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of field galaxies (mostly blue/late type) is higher in the past, almost any of the mechanisms 
mentioned above can explain our observational results. Furthermore, although we discussed 
about cluster specific phenomena, it is also known that galaxies in the field region evolve as 
a function of redshift as well. (e.g. Hammmer et al. 1996; Lilly et al. 1996; Balogh et al. 
1997, 2002). The evolution of field galaxies needs to be compared with that of cluster galaxies 
further in detail. Therefore, it is still an open question what mechanism causes spectral and 
morphological evolution of cluster galaxies. 

The finding of a 30% change of the fraction during the look back time of ~3.5 Gyr could 
also give us an additional hint in finding an underlying physical process. If the gas in spiral 
galaxies is removed very efficiently by some physical processes (e.g. ram-pressure stripping) or 
consumed rapidly by star formation, the spiral arms will disappear after several disk rotation 
periods, ~ 1 Gyr (Sellwood & Carlberg 1984). Interaction/merger processes are quicker than 
gas removal processes (~0.5 Gyr; Mihos 1995). Moore et al.'s (1996) simulation showed that the 
galaxy harassment phase lasts for several Gyr. Kodama et al. (2001) used the phenomenological 
simulations to show that the timescale of the morphological transformation from spiral to SO 
is 1~3 Gyr. For spectral change, Shioya et al. (2002) showed that a disk needs 2-3 Gyr after 
the removal of gas (or truncation of star formation) to show a k spectrum. Poggianti et al. 
(1999) compared the spectral and morphological properties of cluster galaxies and suggested 
that the timescale of the morphological transition is longer than that of the spectral transition. 
This difference in timescale is interesting since if one process is significantly quicker than the 
other, we might be able to see the time difference in the decreases in the fraction of late type 
galaxies and blue galaxies, which will provide a strong constraint in the evolution history of 
the Butcher-Oemler galaxies. In Figure 1, we see a ~30% of change in both the photometric 
and morphological Butcher-Oemler effect between z=0.02 and z=0.30 (~3.5 Gyr). The scatter 
in our measurement, however, is considerable and our choices of criteria between late and early 
type galaxies do not necessarily coincide with each other. It is thus not straightforward to 
convert the information to the time scale of the responsible physical process. In addition, to 
understand change in fraction of morphological and spectral late-type galaxies, the change in 
infalling rate of field galaxies needs to be understood as well. Since computer simulations have 
recently made dramatic progress, in the near future it will become possible for state-of-the-art 
simulations to simulate both dynamical and spectral evolutions of cluster galaxies, plus infalling 
rate of field galaxies in order to compare the results with the observed trend. For example, such 
a simulation can be done by combining dynamical simulations (e.g. Evrard 1991, Kauffmann 
et al. 1995, Bekki, Shioya & Couch 2001, Vollmer et al. 2001, Bekki et al. 2002) with cluster 
phenomenological simulations (e.g. Abraham et al. 1996, Fujita 1998,2001, Balogh et al. 1999, 
Stevens, Acreman, & Ponman 1999, Balogh, Navarro, & Morris 2000, Kodama & Bower 2001). 
Figure 1 in this work provides the interesting observational data to tackle with using such a 
simulation of cluster galaxy formation. 
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4-7. Richness Dependence 

In Section 3.6, we observe the tendency of richer clusters to have smaller fractions of 
late type galaxies, by measuring the residuals from the best fit relations as a function of cluster 
richness. Our result is consistent with Margoniner et al. (2001), who used a similar optical 
richness to find that poorer clusters tend to have larger blue fractions than richer clusters at the 
same redshift. Figure 4, however, still shows a significant amount of scatter, which might be 
suggesting the existence of another factor in determining the blue fraction in addition to redshift 
and richness. The dependence of the late type fraction on cluster richness, however, provides 
another hint on the underlying physical processes. Since ram pressure is stronger in clusters with 
higher temperature at the same gas density, Fujita & Nagashima(1999) theoretically predicted 
that if ram pressure is the only mechanism responsible for the evolution of galaxies in clusters, 
the fraction of blue galaxies will always be higher in lower X-ray luminosity clusters, which 
usually have low temperatures. Our data shown in Figure 4 are consistent with the prediction 
from their ram pressure stripping model. Although our richness is from numbers of galaxies in 
optical imaging data, it is reasonable to assume it correlates well with X-ray luminosity (Bahcall 
1977; Bower et al. 1994). Then, the optical richness can be related to the gas temperature 
through the well known Lx — T relation (Mitchell, Ives, and Culhane 1977; Henry & Arnaud 
1991; Edge & Stewart 1991; David et al. 1993; White, Jones, and Forman 1997; Allen & Fabian 
1998; Markevitch 1998; Arnaud & Evrard 1999; Jones & Forman 1999; Reichart, Castander, 
& Nichol 1999; Wu, Xue, and Fang 1999; Xue & Wu 2000; and see the references therein). 
In a simple estimation, ram pressure is proportional to pv 2 . Lx is proportional to p 2 . From 
the virial theorem, v 2 oc T. The Lx — T relation studied by Xue & Wu (2000) is Lx oc T 2 8 . 
Therefore, ram pressure is proportional to ~ L x m - Combined with an assumption that optical 
richness scales with X-ray luminosity (see, e.g. Bahcall 1974, Jones & Forman 1978, Bower et 
al. 1994, and Miller et al. in preparation), Figure 4 provides another hint that ram pressure 
stripping induces the evolution of cluster galaxies. 

In the literature, however, the dependence of blue fractions on cluster richness has been 
controversial. Bahcall (1977) studied 14 X-ray clusters and found that the fraction of spiral 
galaxies decreases with increasing X-ray luminosity. Lea & Henry(1988) observed 14 clusters in 
X-ray and found that the percentage of blue objects in the clusters seems to increase with the 
X-ray luminosity. On the other hand, Fairley et al. (2002) studied eight X-ray selected clusters 
and did not find any dependence of blue fractions on X-ray luminosities. Balogh et al. (2002) 
studied 10 clusters at z=0.25 with low X-ray luminosity and found similar morphological and 
spectral properties of galaxies compared with clusters with high X-ray luminosity (Balogh et al. 
1997). In all cases, the results were based on a small sample of clusters. We also point out that 
although our results are consistent with a ram-pressure stripping model, there is a possibility 
that the other mechanisms could explain the phenomena. For example, richer clusters might 
have higher rate of merger/interaction due to their higher galaxy density. The same argument 
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holds true for galaxy harassment. Thus, more study is needed to conclude about the physical 
mechanism causing the phenomena. In the near future, confirming the richness dependence 
using X-ray luminosities or velocity dispersions with a larger sample of clusters would offer us 
more insight on the subject. 

5. Conclusions 

In this paper, we have investigated the fraction of late type galaxies in four different 
ways using one of the largest, most uniform samples of 514 clusters between 0.02<z<0.3 from 
the SDSS Cut & Enhance galaxy cluster catalog. All the clusters selected here have more 
than 25 member galaxies within 0.7 Mpc and between M r = —24.0 and M r = —19.44 after 
statistical local fore/background subtraction. The following four different ways to estimate the 
fractions of late type galaxies are adopted: restframe g — r color (a classical Butcher-Oemler 
estimator), u — r color, concentration index and de Vaucouleur/exponential profile fit. The 
last three parameters are known to correlate well with galaxy morphologies (Shimasaku et al. 
2001, Strateva et al. 2001). In all four cases, we observe an increase of the fraction of late type 
galaxies toward higher redshift with a significance of more than 99.99% (Table 1). We draw 
the following conclusion from this work. 

1. We confirm the presence of the Butcher-Oemler effect using g — r color. The Butcher- 
Oemler effect is real and exists in our sample clusters as seen in the lower left panel of Figure 
1. The slope of the increase is consistent with previous work although the scatter in the blue 
fraction is considerable. Previous work also noted a large scatter in the fraction of blue galaxies. 
The fraction of late type galaxies also shows a similar increase when we use a, u — r color cut. 

2. We observe a morphological Butcher-Oemler effect as an increase of late type galaxies 
toward higher redshift, using pure morphological parameters such as a concentration parameter 
and de Vaucouleur/exponential profile fit. The rates of increase are consistent with previous 
work on the spiral to SO transition, albeit with considerable scatter (Figure 1). The increase 
is also in agreement with the original Butcher-Oemler effect from g — r color. Our results are 
consistent with the evolutionary scenario proposed by Dressier et al. (1997), Smail et al. (1997), 
Couch et al. (1998), and Kodama & Smail (2002), in which there is a progressive morphological 
conversion in clusters from spirals into E/SO's. 

3. We find a slight tendency for richer clusters to have lower values of the late type 
fraction (Figure 4). This trend agrees with the ram pressure stripping model proposed by 
Bahcall (1977) and Fujita et al. (1999), in which galaxies in richer clusters are more affected 
by ram pressure due to their high temperature. 

Although our results 1,2, and 3 are all consistent with a ram-pressure stripping model, 
there still remains a possibility that other physical mechanisms are responsible for the evolution 
of cluster galaxies. Thus, further study is needed both theoretically and observationally to reveal 
the underlying physical mechanism responsible for the evolution of cluster galaxies. Since this 
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work is based on only 5% of the whole SDSS data, an increase in the data will improve the 
statistical accuracy as the SDSS proceeds. Extending the work to higher redshifts using 4-8 m 
class telescopes will offer more insight on the origin and evolution of cluster galaxies. 
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Appendix 1. Varying Radius 

In Section 3, we used a fixed 0.7 Mpc radius to measure blue/spiral fractions among 
cluster galaxies since it was difficult to measure virial radius for relatively poor clusters in our 
sample from the SDSS imaging data. However it is known that virial radius changes according 
to cluster richness; i.e. richer clusters have larger virial radius than poorer clusters. Therefore 
using a fixed radius could bring some bias associated with cluster richness. In this section 
we try to rectify this problem using cluster richness to calculate virial radius under a simple 
assumption. We assume that our cluster richness (number of galaxies between M r * = —24.0 
and M r = —19.44 within 0.7 Mpc after fore/background subtraction) is proportional to volume 
of a cluster, and therefore proportional to radius 3 . Since richness is a relatively easy parameter 
to measure from the imaging data, we use the following equation to calculate radius for each 
cluster. 

radius = 0.7 x (Richness /32) X/ ' 3 (Al) 

where median richness of our sample cluster is 32. The coefficient of the equation is 
adjusted so that median clusters in our sample have radius of 0.7 Mpc, which corresponds to 
the mean radius used in Butcher et al. (1978, 1984) and Margoniner et al. (2000, 2001). The 
distribution of radius calculated in this way is presented in figure 12. As expected it has a peak 
at 0.7 Mpc. Using this varying radius, we re-calculated all figures in Section 3. Sample clusters 
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are still required to have more than 25 galaxies after fore/background subtraction within the 
new radius. Therefore the number of sample clusters are somewhat reduced to 413 clusters. 
Results are presented in figures 13-15. Reassuringly, all figures have the same trend as presented 
in Section 3. Therefore the discussion in Section 4 still holds. Although it is ideal to use virial 
radius to measure blue/spiral fractions of clusters, we regard that our analysis using fixed 0.7 
Mpc radius is not hampered to the extend where our main conclusions change. 

Appendix 2. Errors on Blue/Late Type Fractions 

In this section we summarize how we derived eq. (2) to calculate errors on blue/late 
type fractions. As a starting point, we assume the following. 

• Number of galaxies in a cluster region (N^j) follows Poisson statistics. 

• Number of galaxies in a certain area of field region (Nj U ) follows Poisson statistics. 

• A* + ^- and Nf are independent of each other. 

• Number of blue/late type galaxies in a cluster region (N l c +f) is strongly correlated with 
number of all galaxies in that region (N^A. 

• Number of blue/late type galaxies in a certain field region (N l f ate ) is strongly correlated 
with number of all galaxies in that region (Nf 1 ). 

And we clarify the definition of our notation. In this appendix, 5 A means a deviation 
of a sampled value from an expectation value, E(A). 
5A=A — E(A), where A is each data value. 
E(A) and 5 A satisfy the following relations. 

E(5A)—0, E(5A) 2 =a 2 , and if A and B are independent, E(5A5B)=0. Note that the equation 
2 is not a deviation of a single sample but the expectation value estimated from the sample, 
and should be written as E(5f 2 ate ) if we write rigidly. 

Under these assumptions, the error of late type fraction, Sfi ate , become 

SfLe = (X/Y) 2 x ((5X 2 /X 2 ) + (5Y 2 /Y 2 ) - 2(SXSY/XY)), (A2) 

where X=N l c fj - Nf te , Y=Ng f - Nf, and f late = X/Y. 

Since Nf.f and Nf follow Poisson statistics, and they are independent. 

E{5Y 2 ) = Ng f + Nf (A3) 

Similarly, when N£+f follows Poisson statistics, N™jH also follows Poisson statistics since N l c +J- ~ 
Ng f x f late . Therefore, 

E{5X 2 ) = N l £ e f + N l f ate . (A4) 

Deriving the cross term at the end of the equation is not so straightforward. The cross 
term is expanded as 

5X5Y = 8(N l c f e f - N l f e )8(N^ f - Nf) (A5) 
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= (5(N l c f f ) - 5(N l f ate ))(5(N^ f ) - 5(Nf 1 )) (A6) 

= 8(N l c f e f )8(N^ f ) - 5(N l f ate )5(Nf_l f ) - 5(N l c a ^ e f )5(Nf) + 5(N l f ate )5(Nf l ) (A7) 
Since we assume N l j ate and N^j, N£+f and N l j ate are both independent, 

E(-5(N l / te )5(N^ f )) = (A8) 

and 

E(-8(N l c f e f )8(Nf)) = 0. (A9) 
Therefore, 

E{5X5Y) = E{5{N l c a ^ e f )5{N^ f ) + 5(N l f e )5(Nf)) (A10) 
Since we assume that N 1 ^ and N^+t, or N l ^ ate and NJ U strongly correlate, we can 



approximate that 



E{6{N%})6(N£ f )) = tr{N%fM*£f) = V^V^ ( AU ) 

and, 



E{S(N l f ate )5(Nf)) = a(N l f ate )a(Nf) = ^Nf te ^Nf (A12) 
Therefore we obtain, 

E(5X5Y) = y/N l £}y/Ng f + ^Nf te ^Nf ( A13) 

By substituting eq. (A13) for 5X5Y in eq. (A2) , we derive eq. (2). 

However this is not the only way to estimate the error. Actually, the correlation between 
jyiate anc j jyail j g no ^. go obvious since late type fraction in the field and that in the cluster 
region might be different. Although we regard the difference is so small that we can assume 
the eq.(Al2), if we assume that N l f ate , Nf rly (= Nf-N l f ate ), N l c f^, and N™ r j y (= N^ f -N l c f e f ) 
are independent, we derive, 

E(6(N%})6{Ng f )) = E(5(N l c ^ f )(5(N l c f f ) + 8(N e j} y ))) = a(A^) 2 = N l c f f (AU) 

and, 



S(N l f ate )S(Nf) = a(N l f ate ) 2 = N\ ate . 
Then, the expectation value of 5X5Y becomes, 

E(5X5Y) 
In this case, eq. (2) becomes, 



i\rlate i i\rlate 



(A15) 



(A16) 



E(Sfi 



late j 



: flate X 



N late + N late 



late 
c+f 



jylate- 



+ 



Ng f + Nf 



2{Nl ate 



c+f 



Nf te ) 



) 2 (N™ f -Nf) 2 (N l c f e f -Nf te )(N^: f -Nf) 



-(A17) 
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z 



Fig. 1. Photometric and morphological Butchcr-Ocmler effect from the 514 SDSS Cut & Enhance clus- 
ters, /j, fcini fexp and f u - r are plotted against redshift. The dashed lines show the weighted least-squares 
fit to the data. The stars and solid lines show the median values. The median values of errors are shown 
in the upper left corners of each panel. The Spearman's correlation coefficients are shown in Table 1. 
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Blue/late type fraction 



Fig. 2. Normalized distributions of late type fractions (fb, fan, fex P and f u -r)- Dashed lines show distri- 
butions of lower redshift clusters (z<0.15) and solid lines show ones of higher redshift clusters (0.15<z<0.3). 
The results of Kolomogorov-Smirnov tests are shown in Table 2. In all cases, Kolomogorov-Smirnov tests 
show the two distributions are significantly different. 
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Fig. 3. Richness distribution as a function of redshift. Richnesses are measure as the number of galaxies 
brighter than M r = —19.44 within 0.7 Mpc after fore/background subtraction. 
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Fig. 4. The difference of the late type fractions from the best fit lines as a function of redshift are plotted 
against cluster richnesses. Solid lines and stars show the median values. 
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Elongation 

Fig. 5. The difference of late type fractions from the best fit lines as a function of redshift are plotted 
against cluster elongation, which was measured as a ratio of major axis to minor axis on a cluster detection 
map of Goto et al. (2002a). Solid lines and stars show median values. 
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Fig. 6. The distribution of the inverse of the concentration index (Cm), defined as the ratio of Petrosian 
50% flux radius to Petrosian 90% flux radius. The solid line shows the distribution of Cj„ measured in 
the g band image. The dashed line shows the distribution of C;„ measured in the r band image. The 
difference between the g band and r band is marginal, assuring our usage of r band Ci n in the right upper 
panel of Figure f from z=0.02 to z=0.3. The statistics are summarized in Table 4. 
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Fig. 7. The dependence of fa n on seeing. Open squares and solid lines show the distribution and 
medians of low z clusters (z<0.15). Filled triangles and dashed lines show the distribution and medians of 
high z clusters (0.24<z<0.3). The median bins are chosen so that equal numbers of galaxies are included 
in each bin. 
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Fig. 8. The seeing distribution of all galaxies brighter than r=21.5. 87% of all galaxies have seeing better 
than 2.0 arcsec. 
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Fig. 9. The dependence of / exp on seeing. Open squares and solid lines show the distribution and medians 
of low z clusters (z<0.15). Filled triangles and dashed lines show the distribution and medians of high z 
clusters (0.24<z<0.3). Median bins are chosen so that equal numbers of galaxies are included in each bin. 
1 a errors shown as vertical bars are more dominant. There is no significant trend with seeing. 
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Late type fraction 



Fig. 10. Various systematic tests. Solid lines show distributions for 2.1-2.21 Mpc annular 
fore/background subtraction. Dashed lines show distributions for a global background subtraction. Dotted 
lines show distributions using 0.7/(l+z) Mpc radius assuming a standard cold dark matter cosmology. 
Long dashed lines show distributions using the brightest galaxy position as a cluster center. In none of the 
cases does a Kolomogorov-Smirnov test show significant difference between the distributions (significance 
to be different is less than 26% in all cases). 
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Fig. 11. Comparison with late- type fraction from spectroscopy. Late-type fractions measured using 
spectroscopic data are plotted against that from imaging data for three clusters with z<0.06 (ABELL 295, 
RXC J01 14.9+0024, and ABELL 957). Dashed lines are drawn to guide eyes. All points agree with each 
other within the error. 
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Fig. 12. Distribution of varying radius to measure blue/spiral fractions. It has a peak at 0.7 Mpc. 
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Fig. 13. The same as figure 1, but measured with varying radius. 
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Fig. 14. The same as figure 4, but measured with varying radius. 
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Fig. 15. The same as figure 5, but measured with varying radius. 
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Table 1. Spearman's correlation coefficients between z and fractions of late type galaxeis. 514 clusters with richness>25 
are chosen as a sample. 





Correlation coefficient 


Significance 


N clusters 


fb 


0.238 


4.4xl0" 8 


514 


fu — r 


0.234 


7.6xl0~ 8 


514 


fexp 


0.194 


9.6xl0~ 6 


514 


fdn 


0.223 


2.9xl0~ 7 


514 



Table 2. Significances in Kolomogorov-Smirnov tests between distributions for z<0.15 and 0.15<z<0.3. In all cases, 
Kolomogorov-Smirnov tests show the distributions for the lower redshift sample and the higher redshift sample arc signifi- 
cantly different. 

Significance 

f b 2.9xl0- 3 

f u _ r l.OxlO" 3 

f exp 3.4xl0- 4 

fan 2.9X10- 3 



Table 3. Scatters in late- type fractions around the best-fit line are compared with median errors of late- type fraction 
calculated with equation (2) 





fb 


fu—r 


fexp 


fdn 


Real scatter (la) 


0.169 


0.183 


0.171 


0.163 


Error estimate 


0.078 


0.069 


0.050 


0.089 



13 



Table 4. Change in the fraction of galaxies with d„ >0.4 (late type) in two different Glters(g,r). 



band 


N(C m >0.4) 


N(total) 


Percentage(%) 


9 


802 


1336 


60.0 


r 


787 


1336 


58.9 


Difference 


15 


1336 


1.1 



Table 5. Change in the fraction of galaxies with exponential fit likelihood greater than dc Vaucoulcur likelihood (late type) 
in two different filters(g, r). Since we discard the galaxies with the same likelihood in this analysis, the total number of 
galaxies in the sample are different in g and r. 



band 


N(late) 


N(late 4 


- early) Percentage(%) 


9 


503 


804 


62.6 


r 


476 


792 


60.1 


Difference 






2.5 
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